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0. Canned Diatribe, Introduction, and Apologies 

Following the discovery of spacetime anomaly cancellation in 1984 string theory has 
undergone rapid development in several directions. The early hope of making direct con- 
tact with conventional particle physics phenomenology has however long since dissipated, 
and there is as yet no experimental program for finding even indirect manifestations of 
underlying string degrees of freedom in nature @. The question of whether string theory 
is "correct" in the physical sense thus remains impossible to answer for the foreseeable 
future. Particle/string theorists nonetheless continue to be tantalized by the richness of 
the theory and by its natural ability to provide a consistent microscopic underpinning for 
both gauge theory and gravity. 

A prime obstacle to our understanding of string theory has been an inability to pene- 
trate beyond its perturbative expansion. Our understanding of gauge theory is enormously 
enhanced by having a fundamental formulation based on the principle of local gauge in- 
variance from which the perturbative expansion can be derived. Symmetry breaking and 
nonperturbative effects such as instantons admit a clean and intuitive presentation. In 
string theory, our lack of a fundamental formulation is compounded by our ignorance of 
the true ground state of the theory. Roughly two years ago, there was some progress 
towards extracting such nonperturbative information from string theory, at least in some 
simple contexts. The aim of these lectures is to provide the conceptual background for 
this work, and to describe some of its immediate consequences. 
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In string theory we wish to perform an integral over two dimensional geometries and 
a sum over two dimensional topologies, 



where the spacetime physics (in the case of the bosonic string) resides in the conformally 
invariant action 



Here fj,, v run from 1, . . . , D where D is the number of spacetime dimensions, G^ V {X) is the 
spacetime metric, and the integral T>g is over worldsheet metrics. Typically we "gauge-fix" 
the worldsheet metric to g ab = e v 5ab, where f is known as the Liouville field. Following 
the formulation of string theory in this form (and in particular following the appearance 
of ||), there was much work to develop the quantum Liouville theory (some of which is 
reviewed in section 4 here), and conformal field theory itself has been characterized as 
"an unsuccessful attempt to solve the Liouville theory" @. Evaluation of the partition 
function Z above without taking into account the integral over geometry, however, does 
not solve the problem of interest, and moreover does not provide a systematic basis for a 
perturbation series in any known parameter. 

The basic idea of |3|-[5[ relied on a discretization of the string worldsheet to provide a 
method of taking the continuum limit which incorporated simultaneously the contribution 
of 2d surfaces with any number of handles. At one fowl swoop, it was thus possible not 
only to integrate over all possible deformations of a given genus surface (the analog of 
the integral over Feynman parameters for a given loop diagram), but also to sum over 
all genus (the analog of the sum over all loop diagrams). This would in principle free 
us from the mathematically fascinating but physically irrelevant problems of calculating 
conformal field theory correlation functions on surfaces of fixed genus with fixed moduli 
(objects which we never knew how to integrate over moduli or sum over genus anyway). 
The progress, however, is limited in the sense that these methods only apply currently for 
non-critical strings embedded in dimensions D < 1 (or critical strings embedded in D < 2), 
and the nonperturbative information even in this restricted context has proven incomplete. 
Due to familiar problems with lattice realizations of supersymmetry and chiral fermions, 
these methods have also resisted extension to the supersymmetric case. 

The developments we shall describe here nonetheless provide at least a half-step in 
the correct direction, if only to organize the perturbative expansion in a most concise 
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way. They have also prompted much useful evolution of related continuum methods. Our 
point of view here is that string theories embedded in D < 1 dimensions provide a simple 
context for testing ideas and methods of calculation. Just as we would encounter much 
difficulty calculating infinite dimensional functional integrals without some prior experience 
with their finite dimensional analogs ||, progress in string theory should be aided by 
experimentation with systems possessing a restricted number of degrees of freedom. 

These notes have been confined in content essentially to the four lectures actually 
given, in order to keep them reasonably short and accessible. (Other review references 
on the same general subject are |||10|)- This means that we stop well short of some of 
the more interesting recent developments in the field (some of which were covered by later 
lecturers at this school), including the application of the critical D = 2 dimensional models 
to address issues of principle such as topology change in 2d quantum gravity, and their 
relation as well to recent work on D = 2 black holes in string theory. We shall present no 
formal conclusions here other than to note that the subject remains in active development, 
and we have tried at various points in the text to draw attention to issues in need of further 
understanding. 

1. Discretized surfaces, matrix models, and the continuum limit 

1.1. Discretized surfaces 

We begin by considering a U D = dimensional string theory", i.e. a pure theory of 
surfaces with no coupling to additional "matter" degrees of freedom on the string world- 
sheet. This is equivalent to the propagation of strings in a non-existent embedding space. 
For partition function we take 



where the sum over topologies is represented by the summation over h, the number of 



general. The most progress in the continuum has been made via the Liouville approach 
which we briefly review in section 4. If we discretize the surface, on the other hand, it 
turns out that ([T[l]) is much easier to calculate, even before removing the finite cutoff. We 
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consider in particular a "random triangulation" of the surface [11], in which the surface is 
constructed from triangles, as in fig. 1. The triangles are designated to be equilateral,0 so 
that there is negative (positive) curvature at vertices i where the number iVj of incident 
triangles is more (less) than six, and zero curvature when iVj = 6. The summation over 
all such random triangulations is thus the discrete analog to the integral J T*g over all 
possible geometries, 

E fa* - E • d- 2 ) 



genus h random 

triangulations 



Fig. 1: A piece of a random triangulation of a surface. Each of the triangular 
faces is dual to a three point vertex of a quantum mechanical matrix model. 



The discrete counterpart to the infinitesimal volume element yj~g is cr^ = iVj/3, so that 
the total area \S\ = Y^i a i J us ^ counts the total number of triangles, each designated to 
have unit area. (The factor of 1/3 in the definition of Oi is because each triangle has three 
vertices and is counted three times.) The discrete counterpart to the Ricci scalar R at 
vertex i is Ri = 27r(6 — iVj)/iVj, so that 

I y/g R ->^24ir(l - Ni/6) = 4tt{V - \F) = 4tt(V - E + F) = An X ■ 
J i 

1 We point out that this constitutes a basic difference from the Regge calculus, in which the 
link lengths are geometric degrees of freedom. Here the geometry is encoded entirely into the 
coordination numbers of the vertices. This restriction of degrees of freedom roughly corresponds 
to fixing a coordinate gauge, hence we integrate only over the gauge-invariant moduli of the 
surfaces. 



4 



Here we have used the simplicial definition which gives the Euler character x i n terms of 
the total number of vertices, edges, and faces V, E, and F of the triangulation (and we 
have used the relation 3F = 2E obeyed by triangulations of surfaces, since each face has 
three edges each of which is shared by two faces). 

In the above, triangles do not play an essential role and may be replaced by any set 
of polygons. General random polygonifications of surfaces with appropriate fine tuning of 
couplings may, as we shall see, have more general critical behavior , but can in particular 
always reproduce the pure gravity behavior of triangulations in the continuum limit. 

1.2. Matrix models 



We now demonstrate how the integral over geometry in (|1.1| ) may be performed in 
its discretized form as a sum over random triangulations. The trick is to use a certain 
matrix integral as a generating functional for random triangulations. The essential idea 



goes back to work [[T2[ on the large iV limit of QCD, followed by work on the saddle point 



approximation |13 



We first recall the (Feynman) diagrammatic expansion of the (O-dimensional) field 
theory integral@ 



dip -ip 2 /2 + A<//4! 
: e 



^7T 



(1.3) 



where cp is an ordinary real number.tj In a formal perturbation series in A, we would need 
to evaluate integrals such as 



A' 



V/2 ftp 



4! 



Up to overall normalization we can write 







2 k 



dJ 2k 



-p 2 /2 + 



d 2k J 2 /2 

e 



j=o 



dJ 2k 



j=o 



;i.5) 



Since -§je J ^ 2 = Je J ^ 2 , applications of d/dJ in the above need to be paired so that 
any factors of J are removed before finally setting J = 0. Therefore if we represent each 
"vertex" \ip A diagrammatically as a point with four emerging lines (see fig. 2b), then ( |1.4| ) 



2 We apologize for this recapitulation of standard Feynman diagram technology, but prefer to 
keep these notes at least marginally accessible to the mathematics community. 

3 The integral is understood to be defined by analytic continuation to negative A. 



5 



simply counts the number of ways to group such objects in pairs. Diagrammatically we 
represent the possible pairings by connecting lines between paired vertices. The connecting 
line is known as the propagator (iff) (see fig. 2a) and the diagrammatic rule we have 
described for connecting vertices in pairs is known in field theory as the Wick expansion. 



(a) (b) 

Fig. 2: (a) the scalar propagator, (b) the scalar four-point vertex. 

When the number of vertices n becomes large, the allowed diagrams begin to form 
a mesh reminiscent of a 2-dimensional surface. Such diagrams do not yet have enough 
structure to specify a Riemann surface. The additional structure is given by widening the 
propagators to ribbons (to give so-called "thick" graphs). From the standpoint of (|1.3|) , 
the required extra structure is given by replacing the scalar ip by an iV x N hermitian 
matrix M l j. The analog of ( |1.5| ) is given by adding indices and traces: 



-trM 2 /2 d d -trM 2 /2 + tr JM 



e ' M 11 • • • M 



Jl ±V± Jr 



J=0 



8 d trJ 2 /2 

e 



1.6) 



where the source J l j is as well now a matrix. The measure in ( |1.6|) is the invariant AM — 
EL dM\ ]7i<j dReM%- dlmM^ , and the normalization is such that / M e _trM2 / 2 = 1. To 
calculate a quantity such as 



A " t e -* M *i\ tM y , (i.7) 

M 



we again lay down n vertices (now of the type depicted in fig. 3b), and connect the legs 
with propagators (M" l jM k i) = 8j 5j (fig. 3a). The presence of upper and lower matrix 
indices is represented in fig. 3 by the double lines@ and it is understood that the sense of 
the arrows is to be preserved when linking together vertices. The resulting diagrams are 
similar to those of the scalar theory, except that each external line has an associated index 



This is the same notation employed in the large N expansion of QCD [|12[. 
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i, and each internal closed line corresponds to a summation over an index j = 1, . . . , N. 
The "thickened" structure is now sufficient to associate a Riemann surface to each diagram, 
because the closed internal loops uniquely specify locations and orientations of faces. 



(a) 



(b) 



Fig. 3: (a) the hermitian matrix propagator, (b) the hermitian matrix four-point vertex. 

To make contact with the random triangulations discussed earlier, we consider the 
diagrammatic expansion of the matrix integral 



(with M an iV x N hermitian matrix, and the integral again understood to be defined 
by analytic continuation in the coupling g.) The term of order g n in a power series 
expansion counts the number of diagrams constructed with n 3-point vertices. The dual 
to such a diagram (in which each face, edge, and vertex is associated respectively to a dual 
vertex, edge, and face) is identically a random triangulation inscribed on some orientable 
Riemann surface (fig. 1). We see that the matrix integral ( |1.8| ) automatically generates all 
such random triangulationsJl Since each triangle has unit area, the area of the surface is 
just n. We can thus make formal identification with ( |1 . 1|) by setting g = . Actually the 
matrix integral generates both connected and disconnected surfaces, so we have written 
e z on the left hand side of ( |1.8|) . As familiar from field theory, the exponential of the 
connected diagrams generates all diagrams, so Z as defined above represents contributions 
only from connected surfaces. We see that the free energy from the matrix model point of 
view is actually the partition function Z from the 2d gravity point of view. 

5 Had we used real symmetric matrices rather than the hermitian matrices M, the two indices 
would be indistinguishable and there would be no arrows in the propagators and vertices of fig. 3. 
Such orientationless vertices and propagators generate an ensemble of both orientable and non- 
orientable surfaces. 
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There is additional information contained in N, the size of the matrix. If we change 
variables M — > MyN in (|1.8| ), the matrix action becomes iVtr(— |trM 2 + (7trM 3 ), with an 
overall factor of ivJl This normalization makes it easy to count the power of N associated 
to any diagram. Each vertex contributes a factor of N, each propagator (edge) contributes 
a factor of N' 1 (because the propagator is the inverse of the quadratic term), and each 
closed loop (face) contributes a factor of iV due to the associated index summation. Thus 
each diagram has an overall factor 

N V-E+F =N X = N 2-2h ^ (LQ) 

where x is the Euler character of the surface associated to the diagram. We observe that 



the value N = e 7 makes contact with the coupling 7 in (|1.1|) . In conclusion, if we take 
g = e~P and iV = e 7 , we can formally identify the continuum limit of the partition function 
Z in (p..8|) with the Z defined in (|1 . 1| ) . The metric for the discretized formulation is not 
smooth, but one can imagine how an effective metric on larger scales could arise after 
averaging over local irregularities. In the next subsection, we shall see explicitly how this 
works. 

(Actually automatically calculates ( |TTT| ) with the measure factor in ( |1.2| ) cor- 

rected to Y^s \G(S) \ ' wnere 1^(^)1 i s the order of the (discrete) group of symmetries of the 
triangulation S. This is familiar from field theory where diagrams with symmetry result 
in an incomplete cancellation of 1/nVs such as in ( |1.4| ) and (|1 . 7| ) . The symmetry group 
G(S) is the discrete analog of the isometry group of a continuum manifold.) 

The graphical expansion of ( |1.8j ) enumerates graphs as shown in fig. 1, where the 
triangular faces that constitute the random triangulation are dual to the 3-point vertices. 
Had we instead used 4-point vertices as in fig. 3b, then the dual surface would have square 
faces (a "random squarulation" of the surface), and higher point vertices [gk/N k / 2 ~ 1 )txM k 
in the matrix model would result in more general "random polygonulations" of surfaces. 
(The powers of iV associated with the couplings are chosen so that the rescaling M — > 
MyN results in an overall factor of N multiplying the action. The argument leading to 
( |1.9|) thus remains valid, and the power of iV continues to measure the Euler character of 
a surface constructed from arbitrary polygons.) The different possibilities for generating 
vertices constitute additional degrees of freedom that can be realized as the coupling of 2d 
gravity to different varieties of matter in the continuum limit. 



6 Although we could as well rescale M — > M/g to pull out an overall factor of N/g 2 , note that 
N remains distinguished from the coupling g in the model since it enters as well into the traces 
via the N x N size of the matrix. 



8 



1.3. The continuum limit 

From ( |1.9j ), it follows that we may expand Z in powers of N, 

Z(g) = N 2 Z (g) + Z 1 (g) + N~ 2 Z 2 (g) + ... = £ N 2 ~ 2h Z h {g) , (1.10) 

where Zh gives the contribution from surfaces of genus h. In the conventional large iV 
limit, we take N —>■ oo and only Zq, the planar surface (genus zero) contribution, survives. 
Zq itself may be expanded in a perturbation series in the coupling g, and for large order n 
behaves as (see [[14] for a review) 



Z (g) ~ J2 nJ ~ 3 (9/9o) n ~ (9c ~ 9?- 1 ■ (1.11) 



These series thus have the property that they diverge as g approaches some critical coupling 
g c . We can extract the continuum limit of these surfaces by tuning g — > g c . This is because 
the expectation value of the area of a surface is given by 

(A) = (n) = ^- In Z (g) ' 



dg 9- 9 C 

(recall that the area is proportional to the number of vertices n, which appears as the 
power of the coupling in the factor g n associated to each graph). As g — > g c , we see that 
A — > oo so that we may rescale the area of the individual triangles to zero, thus giving a 
continuum surface with finite area. Intuitively, by tuning the coupling to the point where 
the perturbation series diverges the integral becomes dominated by diagrams with infinite 
numbers of vertices, and this is precisely what we need to define continuum surfaces. 

There is no direct proof as yet that this procedure for defining continuum surfaces is 
"correct", i.e. that it coincides with the continuum definition ( |1 . 1| ) . We are able, however, 
to compare properties of the partition function and correlation functions calculated by 
matrix model methods with those (few) properties that can be calculated directly in the 
continuum (for a review, see ||15|1 ). This gives implicit confirmation that the matrix model 
approach is sensible and gives reason to believe other results derivable by matrix model 
techniques (e.g. for higher genus) that are not obtainable at all by continuum methods. 

One of the properties of these models derivable via the continuum Liouville approach 
is a "critical exponent" 7, defined in terms of the area dependence of the partition function 
for surfaces of fixed large area A as 

Z(A) ~ A (7_2)x/2_1 . (1.12) 
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To anticipate some relevant results, we recall that the unitary discrete series of conformal 
field theories is labelled by an integer m > 2 and has central charge D = 1 — 6/m(m+l) (for 
a review, see e.g. jlGy ), where the central charge is normalized such that D = 1 corresponds 
to a single free boson. If we couple conformal field theories with these fractional values 
of D to 2d gravity, the continuum Liouville theory prediction for the exponent 7 is (see 
section 4) 

7 = 1(J3 - 1 - V(£>- !)(£>- 25)) = ~ . (1.13) 

The case m = 2, for example, corresponds to D = and hence 7 = — \ for pure gravity. 
The next case m = 3 corresponds to .D = 1/2, i.e. to a 1/2-boson or fermion. This is 
the conformal field theory of the critical Ising model, and we learn from (|1.13| ) that the 
Ising model coupled to 2d gravity has 7 = — |. Notice that (|1.13j ) ceases to be sensible for 
D > 1. This is the first indication of a "barrier" at D = 1 which will reappear in various 
guises in what follows. 

In section 2 we shall present the solution to the matrix model formulation of the 
problem, and the value of the exponent 7 provides a coarse means of determining which 
specific continuum model results from taking the continuum limit of a particular matrix 
model. Indeed the coincidence of 7 and other scaling exponents (to be defined in sec- 
tion 4) calculated from the two points of view were originally the only evidence that the 
continuum limit of matrix models was a suitable definition for the continuum problem of 
interest. In the past year, the simplicity of matrix model results for correlation functions 
has spurred a rapid evolution of continuum Liouville technology so that as well many 
correlation functions can be computed in both approaches and are found to coincide.B 



7 By way of very superficial overview: following the confirmation that the matrix model ap- 
proach reproduced the scaling results of [17], some 3-point couplings for order parameters at genus 
zero were calculated in from the standpoint of ADE face models on fluctuating lattices. The 



connection to KdV (reviewed in section 3 here) was made in [IS], and then general correlations of 
order parameters (not yet known in the continuum) were calculated in [20|. The first step in the 
calculation of continuum correlators was provided in [21 1, where the free field formulation by zero 



mode integration of the Liouville field was established. This was employed in [22| together with a 
necessary analytic continuation of the scaling parameter to calculate some continuum correlation 
functions: the incorporation of the Liouville mode was shown to cancel the ghastly assemblage 
of T-functions familiar from the conformal field theory result and reproduce the relatively simple 
matrix model result. Additional genus zero correlation functions for D < 1 were then computed in 
[23]. The genus one partition function for the AD series was calculated via KdV methods in pL 
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1.4- The double scaling limit 

Thus far we have discussed the naive N — > oo limit which retains only planar surfaces. 
It turns out that the successive coefficient functions Zh(g) in ( |1.10| ) as well diverge at the 
same critical value of the coupling g = g c (this should not be surprising since the divergence 
of the perturbation series is a local phenomenon and should not depend on global properties 
such as the effective genus of a diagram) . As we shall see in the next section, for the higher 



genus contributions (|1.11|) is generalized to 



z h ( g ) ~ J2 nil ~ 2)x/2 ~ 1 ^/9c) n ~ {g c -g) {2 - 1)x/2 • (1-14) 

n 

We see that the contributions from higher genus (x < 0) are enhanced as g — > g c . This 
suggests that if we take the limits iV — > oo and g — > g c not independently, but together 
in a correlated manner, we may compensate the large N high genus suppression with a 
g — > g c enhancement. This would result in a coherent contribution from all genus surfaces 



To see how this works explicitly, we write the leading singular piece of the Zh(g) as 

Z h (g) ~ f h (g - gJ^M 3 • 

Then in terms of 

K -l = N ( g - g c )V-y)/* , (1.15) 
the expansion (pL . 1 0|) can be rewritten! 

Z = K~ 2 h + fi + K 2 f 2 + ... = J2 K 2h ~ 2 fh ■ (1-16) 

h 

The desired result is thus obtained by taking the limits iV — > oo, g — > g c while holding 



fixed the "renormalized" string coupling k of (|1.15|) . This is known as the "double scaling 
limit" . 



and was confirmed from the continuum Liouville approach in [24]. For D = 1, the matrix model 



approach of |2q,2£| was used in |]27[| (also |2^,^9|) to calculate a variety of correlation functions. 



These were also calculated in the collective field approach [ 30 1 where up to 6-point amplitudes 
were derived, and found to be in agreement with the Liouville results of [23]. 

8 Strictly speaking the first two terms here have additional non-universal pieces that need to 
be subtracted off. 
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2. All genus partition functions 

The large iV limit of the matrix models considered here was originally solved by saddle 
point methods in JT^]. In this section we shall instead present the orthogonal polynomial 
solution to the problem ( ||14|| and references therein) since it extends readily to subleading 
order in N (higher genus corrections). 

2.1. Orthogonal polynomials 

In order to justify the claims made at the end of the previous section, we introduce 
some formalism to solve the matrix models. We begin by rewriting the partition function 
in the form 

Z ^dMe- try(M »=/ndA,A^)e- E - ,/(A - ) , (2.1) 

J i=l 



where we now allow a general polynomial potential V(M). In (|2.1|) , the A^'s are the N 
eigenvalues of the hermitian matrix M, and 

A(A) = n(A, - Xi) (2.2) 

i<j 

is the Vandermonde determinant.i Due to antisymmetry in interchange of any two eigen- 
values, (|2.2|) can be written A(A) = det X{~ x (where the normalization is determined by 
comparing leading terms). In the case N = 3 for example we have 

1 Ai A? 



(A 3 -A 2 )(A 2 -A 1 )(A 3 -A 1 ) = det 1 A 2 A 



1 A 3 \\ 



9 (|2.1| ) may be derived via the usual Fadeev- Popov method: Let Uo be the unitary matrix such 
that M = UqA'Uo, where A' is a diagonal matrix with eigenvalues A^. The right hand side of fl2.1| ) 
follows by substituting the definition 1 = / [J. dAi dU 5{UMU^ - A) A 2 (A) (where J dU = 1). 
We first perform the integration over M, and then U decouples due to the cyclic invariance of the 
trace so the integration over U is trivial, leaving only the integral over the eigenvalues Ai of A. 
To determine A(A), we note that only the infinitesimal neighborhood U = (1 + T)Uq contributes 
to the U integration, so that 

1= / J]dA l dt/5 Jv2 ([/M[/ t -A)A 2 (A)= / dT<5 7V(w - 1) ([T,A'])A 2 (A') . 

i=l ' ! 



Now [T, A']ij = Tij(X'j — X'i), so ( p.2|) follows (up to a sign) since the integration dT above is over 
real and imaginary parts of the off-diagonal Ty's. 
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The now-standard method for solving (^T|) makes use of an infinite set of polynomials 
P n (X), orthogonal with respect to the measure 

/oo 
dA e~ v W P n {\) P m (X) = h n 5 nm . (2.3) 
-oo 

The -P n 's are known as orthogonal polynomials and are functions of a single real variable A. 
Their normalization is given by having leading term P n (X) = A n + . . ., hence the constant 
h n on the r.h.s. of ( |2.3|) . Due to the relation 

A (A) = det A]" 1 = det Pj-x(Xi) (2.4) 

(recall that arbitrary polynomials may be built up by adding linear combinations of 
preceding columns, a procedure that leaves the determinant unchanged), the polynomi- 
als P n can be employed to solve (|2TT|) . We substitute the determinant det Pj-i(Xi) = 
XX— 1)^ Y\ k Pik-i(^k) for each of the A(A)'s in ( |2.1| ) (where the sum is over permutations 
Zfc and (— l) w is the parity of the permutation). The integrals over individual A^'s factorize, 
and due to orthogonality the only contributions are from terms with all Pi(Aj)'s paired. 
There are iV! such terms so (|2.1| ) reduces to 

7 N-l N-l 

e = N\ J] h % = N\ h» J] fk~ k > (2-5) 

i=0 k=l 

where we have defined fk = hk/hk-i- 

In the naive large iV limit (the planar limit), the rescaled index k/N becomes a 
continuous variable £ that runs from to 1, and fk/N becomes a continuous function 
/(£). In this limit, the partition function (up to an irrelevant additive constant) reduces 
to a simple one- dimensional integral: 

k ^° 

To derive the functional form for /(£), we assume for simplicity that the potential 
V(X) in ( |2.3| ) is even. Since the P^s from a complete set of basis vectors in the space 
of polynomials, it is clear that XP n (X) must be expressible as a linear combination of 
lower P^s, XP n (X) = J27=o ai ^M^) (with ai = h^ 1 f e~ v XP n Pi). In fact, the orthogonal 
polynomials satisfy the simple recursion relation, 

XP n = Pn+l + r nPn-l 5 (2-7) 
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with r n a scalar coefficient independent of A. This is because any term proportional to P n 
in the above vanishes due to the assumption that the potential is even, J e~ v X P n P n = 0. 
Terms proportional to Pi for i < n — 1 also vanish since / e~ v P n XPi = (recall XPi is a 
polynomial of order at most % + 1 so is orthogonal to P n for i + 1 < n) . 

By considering the quantity P n XP n -i with A paired alternately with the preceding or 
succeeding polynomial, we derive 



Pn X P n — i T n h n — i h n . 



-V 



This shows that the ratio f n = h n /h n -i for this simple casein is identically the coefficient 
defined by (|2.7|) , f n = r n . Similarly if we pair the A in P^XP n before and afterwards, 
integration by parts gives 

nh n = I e~ v P' XP n = I e~ v P' r n P„._i = r n I e~ v V P n P„_, . (2.1 



n ' ' ~ n I - n 



This is the key relation that will allow us to determine r n . 

2.2. The genus zero partition function 

Our intent now is to find an expression for f n = r n and substitute into ( |2.6|) to 
calculate a partition function. For definiteness, we take as example the potential 



A 4 A 6 
A 3 „, A 5 



(2.9) 



with derivative gV'(X) = A + 2— + 36 ., 

N N z 



The right hand side of ( |2.8| ) involves terms of the form j e~ v A 2p_1 P n P n -\. According to 
( |2.7|) , these may be visualized as "walks" of 2p — 1 steps (p — 1 steps up and p steps down) 
starting at n and ending at n — 1, where each step down from mtom-1 receives a factor 
of r m and each step up receives a factor of unity. The total number of such walks is given 
by (^p 1 ), and each results in a final factor of h n -\ (from the integral J e~ v P n -i P n -i) 
which combines with the r n to cancel the h n on the left hand side of (|2.8|) . For the potential 



(gj), (UD thus gives 



gn = r n + ^r n (r n+1 + r n + r n _i) + ^(10 rrr terms) . (2.10) 



10 In other models, e.g. multimatrix models, /„ = h n /h n -\ has a more complicated dependence 
on recursion coefficients. 
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(The 10 rrr terms start with r n {r\ + r^ +1 + r^ l _ 1 + . . .) and may be found e.g. in [j3Tf .) 

As mentioned before ( |2.6|) , in the large N limit the index n becomes a continuous 
variable £, and we have r n /N — > r(£) and r n ±i/N — > r(£ ±e), where £ = 1/N. To leading 
order in 1/N, (|2.10|) reduces to 



g £ = r + 6r 2 + 306r 3 = W(r) 

2 (2-11) 

In the second line, we have expanded W(r) for r near a critical point r c at which W'\ r = rc = 
(which always exists without any fine tuning of the parameter 6), and g c = W(r c ). We 
see from (|2.11|) that 

r - r c ~ (g c -g£,) 1/2 ■ 

(For a general potential V(A) = f^X]p a p^ 2p i n (|2.9|) , we would have W(r) = 

To make contact with the 2d gravity ideas of the preceding section, let us suppose 
more generally that the leading singular behavior of /(£) (= r(£)) for large N is 

/(O-/c~(0 c -00 -7 (2-12) 

for # near some gr c (and £ near 1). (We shall see that 7 in the above coincides with the 
critical exponent 7 defined in ( |1.12j ).) The behavior of ( |2.6| ) for g near g c is then 



~ y 1 d^ (1 - 0(g c - gO' 1 ~ (i - OCffc - ^)" 7+1 |o + d ^ - 00' 

~(^ c -^)- 7+2 ~E n7 " 3 (^) n - 



-7+1 



(2-13) 

Comparison with (|1.12|) shows that the large area (large n) behavior identifies the expo- 



nent 7 in ( |2.12| ) with the critical exponent defined earlier. We also note that the second 



derivative of Z with respect to x = g c — g has leading singular behavior 

z"~{g a -g)-^~fQ) ■ (2-14) 

From ( |2.12j ) and ( |2.13[ ) we see that the behavior in ( |2.11| ) implies a critical exponent 
7 = —1/2. From ( |1.13| ), we see that this corresponds to the case D = 0, i.e. to pure 
gravity. It is natural that pure gravity should be present for a generic potential. With 
fine tuning of the parameter b in ( |2.9|) , we can achieve a higher order critical point, with 
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W'\ r=rc = W"\ r=rc = 0, and hence the r.h.s. of ( |2.11| ) would instead begin with an (r — r c ) 3 
term. By the same argument starting from ( |2.12| ), this would result in a critical exponent 
7 = —1/3. With a general potential V{M) in ( |2.1| ), we have enough parameters to achieve 
an m th order critical point [[J2] at which the first m—1 derivatives of W{r) vanish at r = r c . 
The behavior is then r — r c ~ (g c — g£,) l t m with associated critical exponent 7 = —1/m. 
As anticipated at the end of subsection 1.2, we see that more general polynomial matrix 
interactions provide the necessary degrees of freedom to result in matter coupled to 2d 
gravity in the continuum limit. 

2.3. The all genus partition function 

We now search for another solution to ( |2.10|) and its generalizations that describes the 
contribution of all genus surfaces to the partition function (|2.6|) . We shall retain higher 



order terms in 1/N in ( 2.10 ) so that e.g. ( 2.11 ) instead reads 



9i = W(r) + 2r(£) (r(£ + e) + r(£ - e) - 2r(0) 

9 (2.15) 
= g c + ±W"\ r=rc {r(0-r c ) + 2r(0(r(£ + e) + r(£ - e) - 2r(£)) + • • • • 

As suggested at the end of section 1, we shall simultaneously let N — > 00 and g — ► g c 
in a particular way. Since g — g c has dimension [length] 2 , it is convenient to introduce a 
parameter a with dimension length and let g — g c = K~ 4 / 5 a 2 , with a — ► 0. Our ansatz 
for a coherent large N limit will be to take e = 1/N = a 5 / 2 so that the quantity k~ 1 = 
(g — g c ) 5 / 4 N remains finite as g — > g c and iV — > 00. 



Moreover since the integral ( |2.6|) is dominated by ^ near 1 in this limit, it is convenient 
to change variables from £ to z, defined by g c — g£ = a z. Our scaling ansatz in this region 



is r(£) = r c + au(z). If we substitute these definitions into fl2.11|) , the leading terms are of 



order a 2 and result in the relation u 2 ~ z. To include the higher derivative terms, we note 
that 

r{£ + e) + r(£ - e) - 2r(0 ~ e 2 -^ = a-^au(z) ~ a 2 W " , 

where we have used e(d/d£) = —ga 1 / 2 (d/dz) (which follows from the above change of 
variables from £ to z). Substituting into ( f2.15| ), the vanishing of the coefficient of a 2 
implies the differential equation 

z = u 2 - ^u" (2.16) 

(after a suitable rescaling of u and z). In ( |2.14j ), we saw that the second derivative of the 
partition function (the "specific heat") has leading singular behavior given by /(£) with 
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£ = 1, and thus by u(z) for z = (g — g c )/a 2 = k -4 / 5 . The solution to (|2.16|) characterizes 
the behavior of the partition function of pure gravity to all orders in the genus expansion. 
(Notice that the leading term is u ~ z 1 ! 2 so after two integrations the leading term in Z 
is z 5 / 2 = k -2 , consistent with ( |1.16p .) 

Eq. (p.!6| ) is known in the mathematical literature as the Painleve I equation. The 
perturbative solution in powers of z~ 5 / 2 = k 2 takes the form u = ^ 1 / 2 (1 — J2k=i u kZ~ 5k ^ 2 ) 1 
where the Uk are all positive.0 This verifies for this model the claims made in eqs. (|1.14|) - 
( |1.16|) of subsection 1.4- For large /c, the Uk go asymptotically as (2/c)!, so the solution 
for u{z) is not Borel summable (for a review of these issues in the context of 2d gravity, 
see e.g. [[33]). Our arguments in section 1 show only that the matrix model results should 
agree with 2d gravity order by order in perturbation theory. How to insure that we are 
studying nonperturbative gravity as opposed to nonperturbative matrix models is still 
an open question. Some of the constraints that the solution to ( |2.16| ) should satisfy are 
reviewed in [31]. In particular it is known that real solutions to ( |2.16[ ) cannot satisfy the 
Schwinger-Dyson (loop) equations for the theory. 

In the case of the next higher multicritical point, with b in ( |2.11|) adjusted so that 
W = W" = at r = r c , we have W(r) ~ g c + ~W"'\ r=rc (r — r c ) 3 + . . . and critical 
exponent 7 = —1/3. In general, we take g — g c = k 2 /^ 7-2 ^ 2 , and e = 1/N = a 2-1 
so that the combination (g — g c ) 1 ~~ / / 2 N = k~ x is fixed in the limit a — > 0. The value 
£ = 1 now corresponds to z = ft 2 '' 7-2 -*, so the string coupling k 2 = z J ~ 2 . The general 
scaling scaling ansatz is r(£) = r c + a~ 2l u(z), and the change of variables from £ to z gives 
e[d/d£) = -ga-"<{d/dz). 

For the case 7 = —1/3, this means in particular that r(£) = r c + a 2 / 3 u(z), k 2 = z -7 / 3 , 



and e(d/d£) = —ga 1 / 3 -^. Substituting into the large iV limit of (|2.10[) gives (again after 



8z 

suitable rescaling of u and z 



u s -uu" -Uu') 2 + au"" , (2.17) 



with a = jq. The solution to ( |2.17| ) takes the form u = 2 1 / 3 (1 + J2k u k z 7k / s ). It turns 
out that the coefficients Uk in the perturbative expansion of the solution to ( |2.17|) are 



11 The first term, i.e. the contribution from the sphere, is dominated by a regular part which 
has opposite sign. This is removed by taking an additional derivative of u, giving a series all of 
whose terms have the same sign — negative in the conventions of ( [2.16| ). The other solution, with 
leading term — z 1 ^ 2 , has an expansion with alternating sign which is presumably Borel summable, 
but not physically relevant. 
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positive definite only for a < jk, so the 3 th order multicritical point does not describe a 
unitary theory of matter coupled to gravity. Although from ( |1.13| ) we see that the critical 
exponent 7 = —1/3 coincides with that predicted for the (unitary) Ising model coupled to 
gravity, it turns out [|3T| , |35|1 that fl2.17p with a = ^ instead describes the conformal field 
theory of the Yang-Lee edge singularity (a critical point obtained by coupling the Ising 
model to a particular value of imaginary magnetic field) coupled to gravity. The specific 
heat of the conventional critical Ising model coupled to gravity turns out (see the next 
section here) to be as well determined by the differential equation ( J2.17| ), but instead with 

For the general m th order critical point of the potential W(r) 1 we have seen that 
the associated model of matter coupled to gravity has critical exponent 7 = — 1/m. With 
scaling ansatz r(£) = r c -\-a 2 t m u(z), we find leading behavior u ~ ^ 1 / m (and Z ~ _ 
k~ 2 as expected). The differential equation that results from substituting the double 
scaling behaviors given before (|2.17[ ) into the generalized version of ( |2.10| ) turns out to be 
the m th member of the KdV hierarchy of differential equations (of which Painleve I results 
for m = 2). In the next section, we shall provide some marginal insight into why this 
structure emerges. 

In the nomenclature of [|36[], so-called "minimal conformal field theories" (those with 
a finite number of primary fields) are specified by a pair of relatively prime integers (p, q). 
(The unitary discrete series is the subset specified by (p, q) = (m + 1, m).) After coupling 
to gravity, these have critical exponent 7 = —2/(p + q — 1). In general, the m th order 
multicritical point of the one-matrix model turns out to describe the (2m — 1,2) model 
(in general non-unitary) coupled to gravity, so its critical exponent 7 = —1/m happens 
to coincide with that of the m th member of the unitary discrete series coupled to gravity. 
The remaining (p, q) models coupled to gravity can be realized in terms of multi-matrix 
models (to be defined in the next section). 



3. KdV equations and other models 

3.1. KdV equations 

We now wish to describe superficially why the KdV hierarchy of differential equations 
plays a role in 2d gravity. To this end it is convenient to switch from the basis of orthogonal 
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polynomials P n employed in the previous section to a basis of orthonormal polynomials 
n n (A) = P n (\) / \fhn~ that satisfy 



oo 

V 



dAe n n n m = S nm . (3.1) 
In terms of the Il n , eq. (|2.7f) becomes 



An n = J -j^- n n + r n y Rn-i = v r n+i n n+ i + v/r^n™-! 

Qnm fT m . 

In matrix notation, we write this as AIT = QH, where the matrix Q has components 

Qnm V m& [ m ,n+l "I - y/^n ^m+l,n ■ (3-2) 

Due to the orthonormality property (|3J]), we see that J e _y AII n II m = Q nrn = Q mn , and 
<5 is a symmetric matrix. In the continuum limit, Q will therefore become a hermitian 
operator. 

To see how this works explicitly |19|J37|| , we substitute the scaling ansatz r(£) = r c + 



a 2 / m u(z) for the m th multicritical model into 

Q^(r c + a 2 /™u(z)) 1 / 2 e * + e 9C (r c + a 2 / m u(z)) 1 / 2 . 
With the substitution e-j^ — > —ga 1 / 771 -^, we find the leading terms 

2/m 

Q = 2r c V 2 + — ■ (u + r c K 2 d 2 ) , (3.3) 

of which the first is a non-universal constant and the second is a hermitian 2 nd order 
differential operator. 

The other matrix that naturally arises is defined by differentiation, 

d 

T^rin = A nm Yl m , (3-4) 

and automatically satisfies [A, Q] = 1. The matrix A does not have any particular symme- 
try or antisymmetry properties so it is convenient to correct it to a matrix P that satisfies 
the same commutator as A. From our definitions, it follows that 

= J -^(U n U m e- v ) A + A T = V'(Q) , 
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where we have differentiated term by term and used J e v X 1 H n H m = (Q^) nm . The matrix 
P = A — TjV' (Q) = \{A — A T ) is therefore anti-symmetric and satisfies 

[p,g] = i. (3.5) 

To determine the order of the differential operator Q in the continuum limit, let 
us assume for example that the potential V is of order 2£, i.e. V = Ylk=o ak ^ 2k • ^ or 
m > n, the integral A mn = J e _v n n ^n m = J e _y Vn n n m may be nonvanishing for 
m — n < 21—1. That means that P mn 7^ for \m — n\ < 2£ — 1, and thus has enough 
parameters to result in a (2£ — l) st order differential operator in the continuum. The single 
condition W = results in P tuned to a 3 rd order operator, and the I — 1 conditions 
W' = ... = W^-V = allow P to be realized (21 — l) st order differential operator. In 
(|3.3|), we see that the universal part of Q after suitable rescaling takes the form Q = d 2 — u. 
For the simple critical point W = 0, the continuum limit of P is the antihermitian operator 
P = d 3 — d}, and the commutator 

l = [P,Q]=4fl£= (5 u a-i«"y (3.6) 

is easily integrated with respect to 2 to give an equation equivalent to ( |2.16| ), the string 
equation for pure gravity. In ( |3.6| ) , the notation R2 is conventional for the first member of 
the ordinary KdV hierarchy. The emergence of the KdV hierarchy in this context is due 
to the natural occurrence of the fundamental commutator relation (|3.5| ), which also occurs 
in the Lax representation of the KdV equations. (The topological gravity approach has as 
well been shown at length to be equivalent to KdV, for a review see [p8| .) 
In general the differential equations 

[P,Q} = 1 (3.7) 
that follow from (|3.5|) may be determined directly in the continuum. Given an operator Q, 



the differential operator P that can satisfy this commutator is constructed as a "fractional 
power" of the operator Q. 

Before showing how this construction works, we first expand slightly the class of 
models from single matrix to multi-matrix models. The free energy of a particular (q — 1)- 



matrix model, generalizing ( |2.1| ), may be written [|39 

Tr 1 -HT t tlVi(M i )-Eti^iMiM i+1 ) 



In / Y[dM, 

J i=i 



In/ J] dA^A(A0e U ' a * [i > + m A(V0 , 



i=l,q-\ 
a=l,N 
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where the Mj (for i = 1, . . . , q — 1) are N x N hermitian matrices, the A^ (a = 1, . . . , N) 
are their eigenvalues, and A(Aj) = n a </?(^! Q ^ — ) ^ s the Vandermonde determinant. 
The result in the second line of (|3.8|) depends on having q's that couple matrices along a 
line (with no closed loops so that the integrations over the relative angular variables in the 
Mi's can be performed.) Via a diagrammatic expansion, the matrix integrals in ( |3.8j ) can 
be interpreted to generate a sum over discretized surfaces, where the different matrices Mj 
represent q — 1 different matter states that can exist at the vertices. The quantity Z in 
( jOf ) thereby admits an interpretation as the partition function of 2d gravity coupled to 
matter. 

Following ||39|| . we can introduce operators Qi and Pi that represent the insertions 
of Xi and d/dA; respectively in the integral ( |3.8| ). These operators necessarily satisfy 
[Pj,Qi] =1. In the iV — ► oo limit, we have seen (following Jl9|]) that P and Q become 
differential operators of finite order, say p, q respectively (where we assume p > q), and 
these continue to satisfy (|3.7| ). In the continuum limit of the matrix problem (i.e. the 
"double" scaling limit, which here means couplings in ( |3lf ) tuned to critical values), Q 
becomes a differential operator of the form 

Q = d q + {v q - 2 {z)A q ~ 2 }+ ■■■ +2v (*), (3-9) 

where d = d/dz. (By a change of basis of the form Q — > f~ 1 (z)Qf(z), the coefficient of 
d 9_1 may always be set to zero.) The continuum scaling limit of the multi-matrix models 
is thus abstracted to the mathematical problem of finding solutions to (|3.7|) . 

The differential equations ( |3.7| ) may be constructed as follows. For p, q relatively 
prime, a p th order differential operator that can satisfy (|3.7| ) is constructed as a fractional 
power of the operator Q of ( |3.9[ ). Formally, a q th root may be represented within an algebra 
of formal pseudo-differential operators (see, e.g. j40[0 as 

oo 

gv^d + ^^d-*} , (3.io) 

i=l 

where d^ 1 is defined to satisfy d -1 / = J^L (— l) J '/w') d -J ' -1 . The differential equations 
describing the (p, q) minimal model coupled to 2d gravity are given by 

[Q+ q , Q] = i, (3.ii) 

where P = Q v J_ q indicates the part of Q v l q with only non-negative powers of d, and is a 
pth orc [ er differential operator. 
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To illustrate the procedure we reproduce now the results for the one-matrix models, 
which can be used to generate (p, q) of the form (21 — 1,2). From (|3.3| ), these models are 
obtained by taking Q to be the hermitian operator 

Q = K = d 2 - u(z) . (3.12) 

The formal expansion of Q^ -1 / 2 = K l ~ 1 / 2 (an anti-hermitian operator) in powers of d is 
given by 

R l-l/2 = d 2l-l _ 2[_1 ^ d 2;_ 3 | + (3 _ 13) 

(where only symmetrized odd powers of d appear in this case). We now decompose 
j^l-l/2 _ j^ 1 ^ i/ 2 _|_^_ 1 / 2 ? w here ifi. 1 ^ 2 = d 2i_1 + • • • contains only non-negative powers 
of d, and the remainder K l _ 1//2 has the expansion 

oo 

K - 1/2 = E {ew-i.d - ^} = {^d- 1 } +0(d" 3 ) + ... . (3.14) 

i=l 

Here we have identified i?z = ei as the first term in the expansion of K l _ . For if 1 / 2 , 

1 12 

for example, we find if + = d and R\ = — u/4. 

The prescription ( |3.11| ) withp = 21 — 1 corresponds here to calculating the commutator 
[K l + , if]. Since if commutes with if' -1 / 2 , we have 

[if+~ 1/2 ,if] = [if,ifi~ 1/2 ] . (3.15) 

But since if begins at d 2 , and since from the l.h.s. above the commutator can have only 
positive powers of d, only the leading (d -1 ) term from the r.h.s. can contribute, which 
results in 

[K l ~ 1/2 , if] = leading piece of [if, 2R t d" 1 ] = AR[ . (3.16) 
After integration, the equation [K l + 1 , if] = 1 thus takes the simple form 

cRi[u]=z, (3.17) 

where the constant c may be fixed by suitable rescaling of z and u (enabled by the property 
that all terms in Ri have fixed grade, namely 21). 

The quantities Ri in (|3. 14 ) are easily seen to satisfy a simple recursion relation. From 



K l+l/2 = KK l~l/2 = K l~l/2 Kj we find 

< +1/2 = \ (K l ~ 1/2 K + KK l ~ 1/2 ) + {R U d} . 
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Commuting both sides with K and using ( |3.16 ), simple algebra gives [IT 



Rl+^jR'/'-uRl-^u'Rt . (3.18) 

While this recursion formula only determines R[, by demanding that the Ri (I ^ 0) 
vanish at u = 0, we obtain 



Rn = - , Rl = U , i?9 = —U 2 u" 

2 ' 4 ' 16 16 

6 32 32 v 2 ; 64 



(3.19) 



We summarize as well the first few K + 



1-1/2 



Kl/ 2 = d 5 - 5 -{u,d 3 } + ^{(3u 2 + u")A} 



(3.20) 



After rescaling, we recognize R3 in ( |3.19| ) as eq. ( |2.17|) with a = i.e. the equation 
for the (2,5) model. In general, the equations determined by ( |3.7| ) for general p, q charac- 
terize the partition function of the (p, q) minimal model (mentioned at the end of section 



2) coupled to gravity. To realize these equations in the continuum limit turns out [fl2| , [43 



to require only a two-matrix model of the type ( |3.8|) . The argument given after (|3.5|) for 



the one-matrix case is easily generalized to the recursion relations for the two-matrix case 
and shows that for high enough order potentials, there are enough couplings to tune the 
matrices P and Q to become p th and q th order differential operators. In subsection 3.2, 
we shall show how to realize a D = 1 theory coupled to gravity in terms of a two-matrix 
model. In |44|], it is argued that one can as well realize a wide variety of D < 1 theories 



by means of a one-matrix model coupled to an external potential. 

3.2. Other models 

As a specific example of a two-matrix model, we consider 

1 dUdV e K N , (3.21) 



where U and V are hermitian N x N matrices and H is a constant. In the diagrammatic 
expansion of the right hand side, we now have two different quartic vertices of the type 
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depicted in fig. 3b, corresponding to insertions of U 4 and V 4 . The propagator is determined 
by the inverse of the quadratic term, 

1 -cV 1 1 (\ c 



-c 1 J 1 - c 2 V c 1 

We see that double lines connecting vertices of the same type (either generated by U 4 or 
V 4 ) receive a factor of 1/(1 — c 2 ), while those connecting U 4 vertices to V 4 vertices receive 
a factor of c/(l — c 2 ). 

This is identically the structure necessary to realize the Ising model on a random 
lattice. Recall that the Ising model is defined to have a spin a = ±1 at each site of 
a lattice, with an interaction Oi<7j between nearest neighbor sites {if)- This interaction 
takes one value for equal spins and another value for unequal spins. Up to an overall 
additive constant to the free energy, the diagrammatic expansion of ( |3.21| ) results in the 
2d partition function 

z= J2 J 3 ^ aiCJj+H ^ ai 

lattices spin 

configurations 

where H is the magnetic field. The weights for equal and unequal neighboring spins are 
e ±/3 , so fixing the ratio e 2 ^ = 1/c relates the parameter c in Q3.21|) to the temperature 
(3. It turns out that the Ising model is much easier to solve summed over random lattices 
than on a regular lattice, and in particular is solvable even in the presence of a magnetic 
field. This is because there is much more symmetry after coupling to gravity, since the 
complicating details of any particular lattice (e.g. square) are effectively integrated out. 

We briefly outline the method for solving ( |3.21|) (see |45|j31| , p5| for more details). By 
methods similar to those used to derive (|2.1] ), we can write fl3.21|) in terms of the eigenvalues 
Xi and yi of U and V, 

z f \ / \ a / \ i r i , -W(xi, yi ) 



Z f 

e = / A(x)A(y) JJ da* dy* e 



where W(x{, yi) = x 2 + yl — 2cxiyi + j^{e H x 4 + e H y 4 )- The polynomials we define for 
this problem are orthogonal with respect to the bilocal measure 



J dxdy e- w ^P n (x) Q m (y) = h n 5 r 
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(where P n ^ Q n for H ^ 0). The result for the partition function is identical to ( pJ^ 

Z 



V X|| //, x I I \ 



and the recursion relations for this case generalize (2~7F) 



V Qmiy^) Qm+1 "I - <?m Qm — 1 ~t~ Qm—3 ■ 

We still have / n = h n /h n -i, and / n can be determined in terms of the above recursion 
coefficients (although the formulae are more complicated than in the one- matrix case). 
After we substitute the scaling ansiitze described in subsection 2.3, the formula for the 
scaling part of / is derived via straightforward algebra. The result is that the specific heat 
u oc Z" is given by (|2.17|) with a = 



Other conventional statistical mechanical models can be formulated on random lattices 
and solved in the continuum limit. The ADE face models (with D < 1), for example, have 
been considered in [J18|. One way of formulating D = 1 is to generalize (p.8|) to an infinite 



line of matrices. In dual form, this is equivalent to strings propagating on a circle of finite 
radius (see e.g. [^5j^6|| ). Another formulation involves letting the index i specifying the 



matrix Mj become a continuous index t G (—00,00). In this limit we trade off matrix 
quantum mechanics for a field theory of matrices theory M(t). This is a problem that 
was originally solved in |T3|], and was used to analyze 2d gravity at genus zero in |Tj| and 



was then applied to higher genus starting in [EREM]. A connection to Liouville theory was 
pointed out in ||7| , and carried further by the free fermion and collective field formulations 



of ^ 

Yet another means of formulating 2d gravity coupled to D = 1 matter is via the 8- 
vertex model, which renormalizes at criticality (the 6- vertex model) onto a single boson at 
finite radius.ll! Since this has not been treated in the literature, we give a quick description 
of the formulation. The simplest vertex models are those for which the degrees of freedom 
are (two-state) arrows that live on links, and are defined on lattices which have four links 
meeting at each vertex. Each possible arrow configuration at a vertex is given a statistical 



12 For an overview geared towards string/particle physicists, see e.g. [48|. On a regular lattice, 
the radius r of the boson (in conventions in which r = l/v2 is the self-dual point) and the 
conventional weights a, b, c of the 6-vertex model are related by cos §7rr 2 = (a 2 + b 2 — c 2 )/2o6. 
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weight, and the partition function is given by summing over all arrow configurations, with 
each assigned an overall weight equal to the product of the statistical weights over the 
vertices. In the 8-vertex model, the vertices are restricted to the set of eight with an even 
number of arrows both incoming and outgoing. In the 6-vertex model, the source and sink 
(all four arrows outgoing or incoming) are excluded, which leaves the four distinct rotated 
versions of fig. 4a, and the two distinct rotated versions of fig. 4b. 



(a) (b) 
Fig. 4: (a) vertex with weight a. (b) vertex with weight c. 

The coupling to gravity is given by summing over random lattices that maintain four 
links at each vertex, but can have arbitrary polygonal faces. It is simple to write down 
a matrix model that generates 6-vertex configurations on random lattices. Rather than a 
hermitian matrix, we employ an arbitrary complex N x N matrix ip = A + iB, where A and 
B are hermitian. The propagator now has an overall orientation, which we identify 

by an arrow on the propagator. (In what follows we suppress the underlying double-lined 
notation of fig. 3.) The graphs of interest are generated by the matrix integral 

e tr(-|^V + a^V 2 + c(y>V)2) ^ 

where the vertices shown in figs. 4a and 4b are assigned weights a and c respectively.li! 

The model has not yet been solved in this formulation except at the analog of the 
Kosterlitz-Thouless point, a = c. At that point we can use the identity 

tr V 2 + V?] =\[{<P + ^) 2 - (¥> - ^) 2 ] 2 = 2 tr(A 2 + S 2 ) 2 



13 On a regular square lattice, the four rotated versions of fig. 4a are further subdivided into 
two mirror reflected pairs, which can be given different weights a and b. On a random lattice 



such distinctions are academic, a property embodied by the cyclicity of the trace in ( |3.22 ), and 
we automatically generate the so-called F-model with a = b. 
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to rewrite the action in terms of the hermitian matrices A, B. By introducing an additional 
integration matrix M, we can reduce the action to terms quadratic in A and B, 

2ctr(A 2 + B 2 ) 2 r tr(-±M 2 + 2^/E(A 2 + B 2 )M) 
e = / e 

Jm 

In this form, the model reduces to a standard transcription of the 0(n) model for n = 2. 
(For general 0(n), A 2 + B 2 is replaced in the above by Yl7=i •) This is reasonable since 
SO (2) is just the circle S 1 normalized to a particular radius. The genus zero solution (due 



to M. Gaudin) is reproduced in [29 



4. Quick tour of Liouville theory 

For completeness, we give here a brief overvie w0 of how the continuum results we 
have used here are calculated. As previously mentioned, the coincidence of these results 
with those of the matrix model approach originally served to give post-facto verification 
of both methods. This section may be considered as an appendix to the preceding three. 

4-1- String susceptibility 7 

We consider the continuum partition function 

VgVX -S M {X;g)-£f&*tJg 

e , (4.1) 



Vol(Diff) 



where Sm is some conformally invariant action for matter fields coupled to a two dimen- 
sional surface £ with metric g, no is a bare cosmo logical constant, and we have symbolically 
divided the measure by the "volume" of the diffeomorphism group (which acts as a local 
symmetry) of £ . For the free bosonic string, we take Sm = §^ / d 2 £ yfg g ab d a X ■ d^X 
where the X(£) specify the embedding of £ into flat .D-dimensional spacetime. 



14 In preparing this section, I may have shamelessly plagiarized some material from a simi- 
lar section in [0] (whose author is consequently responsible for any conceptual errors contained 



herein). Historically, after the work of H some of the results here where derived in [49|, where the 
conformal quantization of Liouville theory was studied (but general correlation functions were not 
calculated) . The quantum Liouville theory was also studied in [|50[ . More recently, the calculation 
of critical exponents in lightcone gauge was carried out in Ji7| (using SL(2, M) current algebra). 
The results were subsequently rederived in conformal gauge in ||51[|, which is the approach we 



follow here since it applies also to higher genus. Reviews of Liouville theory may be found in [52]. 
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To define ( |4.1|) , we need to specify the measures for the integrations over X and g 
(see, e.g. |[53|| ). The measure VX is determined by requiring that JV g 5X e-" 5x ll' = 1, 
where the norm in the gaussian functional integral is given by = J d 2 £ ^fgbX ■ SX. 

Similarly, the measure Vg is determined by normalizing J V g Sg e~?^ S9 ^ 9 = 1, where 
ll^llo = / d 2 £ yfg(g ac g bd + 2g ab g cd ) 8g ab 5g cd , and Sg represents a metric fluctuation at 
some point g- in the space of metrics on a genus h surface. 

The measures T>X and T>g are invariant under the group of diffeomorphisms of the 
surface, but not necessarily under conformal transformations g ab — > e a g ab . Indeed due to 
the metric dependence in the norm ||5X|| 2 , it turns out that 



V e , g X = e^ SLia) V g X , (4.2) 

where 

S L (a) = Jd 2 ^^g ab d a ad b a + Ra + fie^ (4.3) 

is known as the Liouville action. (This result may be derived diagrammatically, via the 
Fujikawa method, or via an index theorem; for a review see ||54j| .) 

The metric measure Vg as well has an anomalous variation under conformal transfor- 
mations. To express it in a form analogous to Q4.2|) , we first need to recall some basic facts 
about the domain of integration. The space of metrics on a compact topological surface E 
modulo diffeomorphisms and Weyl transformations is a finite dimensional compact space 
Aih, known as moduli space. (It is O-dimensional for genus h = 0; 2-dimensional for h = 1; 
and (6/i — 6) -dimensional for h > 2). If for each point r G Aih, we choose a representative 
metric g^, then the orbits generated by the difleomorphism and Weyl groups acting on 
c/ij generate the full space of metrics on E. Thus given the slice g(r), any metric can be 
represented in the form 

r 9 = e^g(r) , 

where /* represents the action of the diffeomorphism / : E — > E. 



Since the integrand of ( ^4.1| ) is diffeomorphism invariant, the functional integral would 
be infinite unless we formally divide out by the volume of orbit of the diffeomorphism 
group. This is accomplished by gauge fixing to the slice g(r); the Jacobian that enters can 
be represented in terms of Fadeev- Popov ghosts, as familiar from the analogous procedure 
in gauge theory. We parametrize an infinitesimal change in the metric as 

s 9 Z z = V 2 ^ , 8g„ = V 2 £z 
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(where for convenience we employ complex coordinates, and recall that the components 
9zz = 9 ZZ are parametrized by e^). The measure T>g at g(r) splits into an integra- 
tion [dr] over moduli, an integration T>ip over the conformal factor, and an integra- 
tion T>£ T>^ over diffeomorphisms. The change of integration variables T>5g zz T>5g zz = 
(det V z det Vg ) £>£ £>£ introduces the Jacobian det V 2 det for the change from 5g to £. 
The determinants in turn can be represented as 

/ d 2 £ ^gb zz V- Z c z -J d 2 £ y/gbzs V 2 c 2 ~ 



detV z detV 5 = / VbVcVbVc e 

(4.4) 

. -5 gh (6,c,6,c) 
P(gh) e 

where T>(gh) = Vb'Dc'Db'Dc is an abbreviation for the measures associated to the ghosts 
b,c,b,c; b zz is a holomorphic quadratic differential, and c 2 (c 2 ) is a holomorphic (anti- 
holomorphic) vector. 

Finally, the ghost measure T>(gh) is not invariant under the conformal transformation 
g — > e a g, instead we have [|6|, |53|j5^1 



2V 9 (gh)=e ^ ^(gh) , (4.5) 



where Sl is again the Liouville action ( |4.3|) . (In units in which the contribution of a single 
scalar field to the conformal anomaly is c = 1, and hence c = 1/2 for a single Majorana- 
Weyl fermion, the conformal anomaly due to a spin j reparametrization ghost is given by 
c = (— 1) F 2(1 + 6j(j — 1)). The contribution from a spin j = 2 reparametrization ghost is 
thus c = —26.) 

We have thus far succeeded to reexpress the partition function (|4.1| ) as 

r -S M - S Rh - ( d 2 £ Jg 

Z = I [dr] V gV > V g (gh) V g X e gh 2w J ^\ 

Choosing a metric slice g = e^g gives 

V eVg ipV eV g{g\Y)V eV gX = j(<P,g) Vg(pV g (gti)T>gX , 

where the Jacobian J((p,g) is easily calculated for the matter and ghost sectors (( |4.2| ) and 
( |4.5|) ) but not for the Liouville mode p>. The functional integral over cp is complicated by 
the implicit metric dependence in the norm 

2 
9 



J d 2 ^(M 2 = / d 2 ^e^ (Sip) 
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since only if the e^ factor were absent above would the T> g ip measure reduce to that of a 
free field. 

In ||51|| , it is simply assumed^ that the overall Jacobian J(ip,g) takes the form of an 
exponential of a local Liouville-like action J d 2 £ \fg (a g ab d a ipdbip + bRip + /ie cv ), where a, 
b, and c are constants that will be determined by requiring overall conformal invariance (c 
is inserted in anticipation of rescaling of ip) . With this assumption, the partition function 
( |4.1|) takes the form 

7 fu it. t» f mt» v ~ S M{X,g) - 5 gh (6,c,6,c; g) 
Z= [dT\Vg(pVg(gh)VgX e 

J (4.6) 

- / d 2 ^ y/g (ag ab d a ipd h ip + bRip + ^) 

■ e 

where the ip measure is now that of a free field. 

The path integral (|4.6| ) was defined to be reparametrization invariant, and should 
depend only on e^g = g (up to diffeomorphism) , not on the specific slice g. Due to diffeo- 
morphism invariance, ( |4.6| ) should thus be invariant under the infinitesimal transformation 

Sg = e{Qg , Sip = -e(£) , (4.7) 



and we can use the known conformal anomalies ( |4.2|) and ( |4.5|) for ip, X, and the ghosts 
to determine the constants a, 6, c. Substituting the variations (JO]) i n (|4.6| ), we find terms 
of the form 

+ ~ b )j d2 ^R £ and (2o-6)y d 2 £y/jjen<p, 

where the D — 26 on the left is the contribution from the matter and ghost measures T> g X 
and Pjj(gh), and the additional 1 comes from the T>§(p measure. Invariance under ( |4.7| ) 
thus determines 



S=?^, a=|S. (4, 



Substituting the values of a, b into the Liouville action in (|4.6|) gives 

h I d2 ^9^^9 ab d a ^pd b ip+ 2 -^-Rip) . (4.9) 



15 Some attempts to justify this assumption may be found in [55]. 
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To obtain a conventionally normalized kinetic term J (<9y?) 2 , we rescale (p — > \J 25-D 
(This normalization leads to the leading short distance expansion <p(z) <p(w) ~ — log (2 — 
w).) In terms of the rescaled (p, we write the Liouville action as 



8tt 

where 



^ / d 2 Zy/jj(g ab da<pd b <p + QR<p} , (4.10) 



(4-11) 



The energy-momentum tensor T = —hd<pd(p + ^-d 2 (p derived from ( [1.10|) has leading short 



distance expansion T(z)T(w) ~ ^CLiouviiie/ {z — w) 4 + . . ., where CLiouviiie = 1 + 3Q 2 . Note 
that if we substitute (|4.11| ) into CLiouviiie an d add an additional c = D — 26 from the matter 
and ghost sectors, we find that the total conformal anomaly vanishes (consistent with the 
required overall conformal invariance). 

It remains to determine the coefficient c in (f4.6j ) . We have since rescaled (p, so we write 
instead e aip and determine a by the requirement that the physical metric be g = ge aip . 
Geometrically, this means that the area of the surface is represented by J (\. 2 ^y/ge aif . a 
is thereby determined by the requirement that e aip behave as a (1,1) conformal field (so 
that the combination d 2 ^e a ^ is conformally invariant). For the energy- momentum tensor 
mentioned after fl4.11| ) , the conformal weightlll of e aip is 

h{e atp ) = h(e a(p ) = -\a{a-Q) . (4.12) 

Requiring that h(e°"^) = h(e aip ) = 1 determines that Q = 2/ a + a. Using ( |4.11| ) and 
solving for a then give 



a 



(V25 — D — VT^D) . (4.13) 



For spacetime embedding dimension d < 1, we find from ( [4 . 1 1| ) and ( [4.13j ) that Q and 
a are both real (with a < Q/2). The D < 1 domain is thus where the Liouville theory is 



Recall that h is given by the leading term in the operator product expansion 
he aip I (z — w) 2 + . . . . Recall also that for a conventional energy-momentum tensor T = —-^dipdip, 
the conformal weight of e xpv is h = h = p 2 /2. 

17 The choice of root for a is determined by making contact with the classical limit of the 
Liouville action. Note that the effective coupling in ( |4.9| ) goes as (25 — D)" 1 so the classical limit 
is given by D — > — 00. In this limit the above choice of root has the classical a — > behavior. 
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well-defined and most easily interpreted. For D > 25, on the other hand, both a and Q are 
imaginary. To define a real physical metric g = e aip g, we need to Wick rotate ip — > —iip. 
(This changes the sign of the kinetic term for (p. Precisely at D = 25 we can interpret 
X° = —iip as a free time coordinate. In other words, for a string naively embedded 
in 25 flat euclidean dimensions, the Liouville mode turns out to provide automatically 
a single timelike dimension, dynamically realizing a string embedded in 26 dimensional 
minkowski spacetime. Finally, in the regime 1 < D < 25, a is complex, and Q is imaginary. 
Sadly, it is not yet known how to make sense of the Liouville approach for the regime of 
most physical interest. We mention as well that so-called non-critical strings (i.e. whose 
conformal anomaly is compensated by a Liouville mode) in D dimensions can always be 
reinterpreted as critical strings in D + 1 dimensions, where the Liouville mode provides 
the additional (interacting) dimension. (The converse, however, is not true since it is not 
always possible to gauge-fix a critical string and artificially disentangle the Liouville mode.) 

It remains to extract the string susceptibility 7 of (|1.13|) in this formalism. We write 
the partition function for fixed A as 

VipVX e~ S^fd 2 ^ v^e QV - A J , (4.14) 

where for convenience we now group the ghost determinant and integration over moduli 
into VX. We define a string susceptibility 7 as in ( |1.14| ) by 

Z(A) ~ A (7 ~ 2)x/2_1 , A -> 00 , 

and determine 7 by a simple scaling argument. (Note that for genus zero, we have Z(A) ~ 
A 7-3 as in ( |1.11| ).) Under the shift <p — > ip + p/a for p constant, the measure in ( [4. 14] ) 
does not change. The change in the action ( fl.lO| ) comes from the term 

Substituting in (|4.14|) and using the Gauss-Bonnet formula j- J d 2 £ y/gR = x together 
with the identity 5(Xx) = S(x)/\X\ gives Z(A) = e~ Qpx / 2a - p Z(e~ p A). We may now 
choose e p = A, which results in 

Z{A) = A~ Qx/2a ~ 1 Z(l) = A (7 ~ 2)x/2_1 Z(l) , 

and we confirm from ( |4.11| ) and ( [4 . 1 3| ) that 

7 = 2 - § = i (D - 1 - y/(D - 25) (D - 1) ) . (4.15) 
a 12 

This result reproduces (|1.13|) , which we used to compare with the result of the matrix 
model calculation (recall that 7 = — 1/m for D = 1 — 6/m(m + 1) ). 
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4-2. Dressed operators / dimensions of fields 

Now we wish to determine the effective dimension of fields after coupling to gravity. 
Suppose that $o is some spinless primary field in a conformal field theory with conformal 
weight ho = h($o) = h($o) before coupling to gravity. The gravitational "dressing" can be 
viewed as a form of wave function renormalization that allows $o to couple to gravity. The 
dressed operator $ = e^$o is required to have dimension (1,1) so that it can be integrated 
over the surface E without breaking conformal invariance. (This is the same argument used 
prior to (|4.13 ) to determine a). Recalling the formula ( f!.12| ) for the conformal weight of 



e^, we find that (3 is determined by the condition 

h - \(3{(3 - Q) = 1 . (4.16) 

We may now associate a critical exponent h to the behavior of the one-point function 
of $ at fixed area A, 

1 f S 

F*{A) = -— J VipVX e" 8^Jd 2 ^e^-A) fd^^e^^-A 1 ^ . (4.17) 

This definition conforms to the standard convention that h < 1 corresponds to a relevant 
operator, h = 1 to a marginal operator, and h > 1 to an irrelevant operator (and in 
particular that relevant operators tend to dominate in the infrared, i.e. large area, limit). 

To determine h, we employ the same scaling argument that led to ( ^.15| ). We shift 
(p — > (p + pjot with e p = A on the right hand side of ( |4.17| ), to find 

A-Qxf2ct-\+P/a 

MA) = A _ Qxl2a _ x F 9 (l) = AV°F 9 (1) , 

where the additional factor of e p ^/ a = A^l a comes from the e^ gravitational dressing of 
$o- The gravitational scaling dimension h defined in fl4.17j ) thus satisfies 

h = l-P/a. (4.18) 

Appealing to the semiclassical argument employed before ( |4.13| ) , we solve ( |4.16| ) for f3 with 
the branch 

P = kQ - \/\Q 2 ~ 2 + 2h o = -^(V2b-D- v/l-^ + 24/io) 
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(for which (3 < Q/2, and (3 — > as D — > — oo). Finally we substitute the above result for 
and the value (p^D for a into flEISD , and findS 



V25-D- y/T=D 



As an example, we apply these results to the minimal models [|36| mentioned at the 
end of section 2. These have a set of operators labelled by two integers p, q (satisfying 
1 < r < q — 1, l<s<p — 1). Coupled to gravity, these operators turn out to have dressed 
conformal weights 

h riS = P + Q ~ lpr ~ qsl l<r<q-l, l<s<p-l, (4.20) 
p + q — 1 

in agreement with the weights determined from the (p, q) formalism discussed in section 3 
for the generalized KdV hierarchy (see e.g. |19| , |56H ) ■ 



More explicitly, we consider the first member of the unitary discrete series, i.e. the D = 
1/2 Ising model, which has (p, q) = (4,3). Before coupling to gravity, critical exponents 
v,ct,/3 can be defined in terms of the divergences of correlation length £ ~ t~ v , specific 
heat C ~ t~ a , and magnetization m ~ t@ with respect to the deviation t = (T — T c )/T c 
from the critical temperature T c . In terms of the conformal weights of the energy and spin 
operators h e and h s , these exponents satisfy v = 2 {i-h ) ; a = ^(1 — ^), P = (2 — a)h s . 
According to ( 4.20 ), the coupling to gravity induces the shifts h e = \ — > |, h s = —> ^, 



which implies corresponding shifts in v, a, and (3. 



We can also substitute = ot{l — 
before ( |4.13| )) to rederive the result h — 
weight" h and the bare weight ho [17]. 



h) from ( |4.18[ ) into ( |4.16| ) and use - \a{a - Q) = 1 (from 
/io = h(l — h)a 2 /2 for the difference between the "dressed 
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